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Abstract 

Dynamical scaling in a two-dimensional lattice model of chaotic maps, in 
contact with a thermal bath at temperature T, is numerically studied. The 
model here proposed is equivalent to a conserved Ising model with couplings 
which fluctuate over the same time scale as spin moves. When couplings 
fluctuations and thermal fluctuations are both important, this model does 
not belong to the class of universality of a Langevin equation known as model 
B; the scaling exponents are continuously varying with T and depend on the 
map used. The universal behavior of model B is recovered when thermal 
fluctuations are dominant. 
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The kinetics of phase separation as a system is quenched below its ordering temperature 
is the subject of much current research In the dynamical scaling regime the ordering 
process is characterized by a single time dependent length R(t) ~ t z representing the average 
domain size 0; the growth exponent z usually does not depend on the dimensionality of the 
system nor on the final temperature of the quenching. In most models with non-conserved 
order parameter z = 1/2, while in scalar models with conserved order parameters z — 1/3 |§. 
A new exponent 9 has been recently introduced by Derrida et al. || with the hope to better 
characterize the process of phase separation: it is related to the persistence probability p(t) 
defined as the probability that the local order parameter at a given point has never changed 
since the initial time; tipically one finds p(t) ~ t~ e Q. 

In a recent paper phase separation mechanisms have been investigated in the framework 
of Coupled Chaotic Lattices ||, reaction-diffusion systems in which chaotic maps are cou- 
pled diffusively ||. In the strong coupling regime these models exhibit nontrivial collective 
behavior, i.e. long-range order emerging out of local chaos. In || the one-body probabil- 
ity distribution functions (pdfs) of local (continuous) variables had two disjoint supports; 
this allowed the unambiguous definition of Ising spin variables. Starting from uncorrelated 
initial conditions, for large coupling values, complete phase ordering was observed; due to 
the deterministic nature of the system, the corresponding Ising spin system was at zero 
temperature. Both the domain size R(t) and the persistence p(t) showed scaling behavior 
and the exponents z and 9 were found to vary continuously with parameters, at odds with 
traditional models. Normal phase ordering behavior, corresponding to the time dependent 
Ginzburg-Landau equation, was recovered in a continuous space limit of the model f7| . These 
findings suggest that the phase ordering properties of multiphase Coupled Chaotic Lattices 
are different from those of most models studied traditionally. It is interesting, therefore, to 
investigate the nonuniversality of scaling exponents in other systems of chaotic maps. 

In this work, we numerically study the dynamical scaling in a lattice model of chaotic 
maps such that the corresponding Ising spin model conserves the order parameter. By 
assuming the system to be in contact with a thermal bath, we also study the competition 

2 



between the autonomous fluctuations of couplings and thermal fluctuations. 

On each site i of a two dimensional square lattice, consider a dynamical system described 
by the variables Xi and yi, which evolve according to the following map: 
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where ^(t) = KJ Xi(t), and cjj(i) = JK(yi(t) — J^x^t)); K, J , and r are free param- 
eters. This map has been studied in [H| in the frame of chaotic neural networks; for a 



suitable choice of the parameters, it has two symmetric chaotic attractors, one with x > 
and the other with x < 0. Hence an Ising spin <Ji(t) = sign[xi(t)} can be associated to 
each dynamical system. At each discrete timestep t we assume the following dynamics, 
consisting of two stages: firstly all the maps are iterated. Secondly, all pairs of nearest 



neighboring maps are sequentially swapped with the following exchange probability [[I I 
Pswap = l/(l + exp{/5A_E t }), where (3 is the inverse temperature, AE t is the change in energy 
that would occur if the maps were exchanged, and the configurational energy is defined as 
E t — — %i(t)xj(t), the sum being over all the nearest neighboring pairs. It is clear that 
this dynamics conserves the spin magnetization by writing E t — — Jij(t)&i(f)&j{t)-i 
with Jij(t) = \xi(t)xj(t)\, we observe that this model is equivalent to a conserved Ising model 
with couplings which fluctuate over the same time scale as spin moves. To study the phase 
ordering process, uncorrelated initial conditions were generated as follows: one half of the 
sites were chosen at random and the corresponding values of x and y were assigned according 
to the invariant distribution of the chaotic attractor with x > 0, while the other sites were 
similarly assigned values with x < 0. Large lattices (up to 1000 x 1000 sites) with periodic 
boundary conditions were used; the persistence p(t) and the domain size R(t) were mea- 
sured. R(t) was estimated by the condition C(R(t),t) = 1/2, where C(r,t) = (a i+r (t)ai(t)) 
is the two point correlation function. The persistence p(t) was measured as the proportion 
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of sites which have not changed a from the initial time. Both R(t) and p(t) were averaged 
over many different samples of initial conditions. 

Let us now discuss our results. In Fig. 1 typical configurations of the ordering system, 
at various times, are shown: the morphology of growing domains is similar to those in the 
conserved Ising model [Q]. At zero temperature {(3 = oo) the system coarsens and exhibits 
scaling behavior. Choosing K = 10, r = 5 and Jo = 0.9, we measure z = 0.11 and 9 = 0.59. 
In Fig. 2a the time evolution of R(t) is shown, while in Fig. 2b the time evolution of the 
persistence p(t) is shown. We remark that the conserved Ising model with spin-exchange 
dynamics (which is recovered by suppressing couplings fluctuations in the present model) 
does not coarsen at zero temperature ||. In the case of finite temperature, we find that 
the growth exponent varies continuously with (3. In Fig. 3 we plot our measures of z for 
some values of the temperature T = 1/(3. We observe that z increases monotonically with 
T and reaches the value 1/3 at T « 0.125; it remains constant until T = 0.32. Therefore, 
in the range T e [0.125,0.32], thermal fluctuations dominate over couplings fluctuations 
and the system belongs to the class of universality of a Langevin equation known as model 
B H, which describes the standard conserved Ising model (when bulk diffusion dominate 
over surface diffusion ||). At high temperature (T > 0.375) the system does not order and 
R(t) tends to a constant after a little transient (z = 0). For temperatures in the range 
[0.32,0.375], our numerical plots of R(t) do not show a neat scaling behavior (filled area 
in Fig. 3): probably, in the thermodynamic limit a critical value (3 C exists, in this range, 
separating the 2=1/3 behavior from the z = one. Further numerical simulations are 
needed to study the phase ordering properties of this system close to (3 C . We note that 
the order of magnitude of (3 C here estimated is consistent with the critical coupling of the 
two-dimensional Ising model (fl c = 0.44) if one approximates the average value of couplings 
as follows: (Jy) ~ (\xi(t)\)(\xj(t)\) = 0.20, the average of x's being calculated over the 
invariant distribution of the chaotic attractor. 

It is worth stressing that simulations of the spin-exchange Ising model, at low tempera- 
ture, do not provide z = 1/3 but a value in the range 0.17 — 0.25, because of corrections due 
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to excess transport in interfaces ||. We cannot exclude that these effects, peculiar to spin- 
exchange dynamics, still play a role, at low temperature, in the present model. However, we 
simulated, at low temperatures, the spin-exchange Ising model and, fitting the domain size 
with the power law form, we measured growth exponents differing from those measured in 
our model (using the same ratio T/T c in the two models). 

Now we consider the scaling function. In the scaling regime the correlation function obeys 
the following scaling form C(r,t) = f(r/R(t)). We remark that recently the universality 
of the scaling function /, with respect to details of the system, has been questioned jl2| . 



In Fig. 4 the scaling collapse of the correlation function is depicted for three values of 
the temperature. At small r/R the scaling function has a linear behavior (Porod's law) 
/ ~ 1 — a r/R; we measure a ~ 0.5, independently of the temperature. On the other hand, 
at larger distances / shows some dependence on the temperature (see Fig. 4): at the moment, 
we have no argument to explain this dependence. Concerning the persistence exponent, at 
finite temperature we measure the block persistence puock{t), which, for a sufficiently large 
block, is much less affected by thermal corrections than p(t) []13|| . We find that also 9 depends 
on the temperature; for example at /3 = 20, using 8x8 blocks, we measure 9 = 1.05. 

We also simulated systems with different maps. We do not report here details and only 
quote the results. We find that, at fixed temperature, the scaling exponents depend on the 
map. Changing the parameters in Eqs. (1-2) and fixing K — 10, r — 5, J Q — 0.95, we 
measure, at zero temperature, z = 0.22 and 9 = 1.57 ; using the one-dimensional map of 
Ref. || (eq. 2), at zero temperature, we find z = 0.07 and 9 = 0.09. In both these two 
cases, however, we verify that the z exponent grows continuously with temperature and, 
when thermal fluctuations are dominant, z is close to 1/3. The scaling function / is found 
to be independent of the map at zero temperature. 

In summary, we have numerically studied the dynamical scaling in a lattice model of 
chaotic maps which is equivalent to a conserved Ising model with fluctuating couplings, in 
contact with a thermal bath at temperature T. When couplings fluctuations and thermal 
fluctuations are both important, this model does not belong to the class of universality 



of model B: the scaling exponents are continuously varying with T and depend on the 
map used. The universal behavior of model B is recovered when thermal fluctuations are 
dominant. Our measures of the scaling function show some dependence on the temperature. 
A similar variation of the scaling exponents with parameters was found in |J , where the Ising 
models associated to the Coupled Map Lattices were not conserving the order-parameter 
and at zero temperature: here we have confirmed this behavior in another class of Chaotic 
Map Lattices, whose corresponding Ising models are conserving the order parameter and in 
presence of thermal fluctuations. 

Finally, we observe that, on the long term, local chaos may be seen as bounded effective 
noise. It follows that a similar behavior might be found in the dynamical scaling of conserved 
coarsening systems in presence of multiplicative noise |14|], i.e. external fluctuations of a 



control parameter in the system. Further investigation, especially at the analytic level, is 
needed to clarify the origin of the nonuniversality here described. 
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Figure Captions 



Figure 1: Snapshots of the ordering system. Black (white) pixels correspond to a — 1 
(-1). Lattice of 100 x 100 sites, (3 — 10 and iterations times t — (a), t — 200 (b), 
t = 25, 000 (c), t = 1, 000, 000 (d). 

Figure 2: Time evolution of the domain size R(t) (a) and persistence p(t) (b) at zero 
temperature. Solid lines are best linear fits. 

Figure 3: The estimated growth exponent z versus the temperature T = 1/(3. 

Figure 4: Scaling collapse of the correlation function for (3 = oo (circles), (3 = 20 (squares) 
and f3 = 10 (triangles). Correlations at eight times, equally spaced in the interval [10 4 , 
4 x 10 4 ], are shown. 
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